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Lanthanum and yttrium compounds with composition RE2M2O7 (RE=Y, La; M = Ti, Zr, Hf) 
have high density and high Z and can be doped with Ce onto the La and Y sites. This makes 
these compounds good candidates for Ce activated scintillator 7-ray detectors particularly for the 
hafnate systems which have a very high density There is disagreement in the literature concerning 
La2Hf 2 07:Ce as it has been reported to show both bright as well as no Ce activated luminescence by 
different experimental groups. We have performed first-principles electronic structure calculations 
of these compounds doped with Ce using the pseudopotential method based on the generalized 
gradient approximation in density functional theory. The positions of the Ce 4/ states relative to 
the valence band maximum and the position of the Ce 5 c? states relative to the conduction band 
minimum (CBM) of the host material are determined. We find, unlike Ce activated La and Y 
compounds where the CBM is typically of La 5c? or Y 4c? character, that, in these systems the CBM 
is predominately of d character on the Ti, Zr, Hf atoms. For all these compounds we also find that 
the Ce 5c? state lies above the CBM which would prevent any luminescence from the Ce site. 



I. INTRODUCTION 

Today Ce 3+ is the most common trivalent lanthanide 
activator for fast and bright scintillation. [T]-[3] Excel- 
lent energy resolution together with high light yield 
and fast decay time has stimulated active interest in 
the development of Ce-doped LaBr 3 for application in 
space missions, [4 medical imaging[5] and radio-isotope 
identification. [6 However, even the efficient rare-earth 
halide scintillators have undesirable attributes like hygro- 
scopic nature and mechanical fragility which pose chal- 
lenges to manufacturing [7H2] and, hence, the quest to 
find new scintillator materials continues. The possibility 
of computational modeling to aid scintillator discovery 
has been the subject of several papers in the past. [T0HT3] 
This is pertinent even more now with the availability of a 
vast literature of known luminescent materials and pow- 
erful computational tools (hardware and algorithms) to 
study underlying physical phenomena and then employ 
them for selecting promising candidate materials. 

We have recently developed a systematic approach 
to select candidate materials for bright Ce scintillation 
based on a set of theoretical criteria derived from stud- 
ies of more than 100 crystalline compounds by first- 
principles electronic structure calculations within the 
framework of density functional theory (DFT).[14 Our 
calculations have been validated for known bright Ce 
scintillators and non-scintillators (compounds not acti- 
vated by Ce 3+ dopant). New Ce-activated scintillator 
materials have also been reported. The objective of our 
theoretical work is to select candidate materials for syn- 
thesis as well as complement experimental work through 
simulations of promising scintillator materials. [15 Since 
Ce 3+ activation has been documented in a variety of host 
materials, [HI [17] the effectiveness of our prescreen selec- 
tion criteria is further crucially tested by how well we are 



able to predict non-scintillators. 

RE 2 M 2 7 (RE = Y, La; M = Ti, Zr, Hf) com- 
pounds have been actively studied for applications in 
nuclear waste removal, [18] catalysis, [19] as hosts for lu- 
minescent centers, [20 and recently in the semiconduc- 
tor industry. [21] Our interest in these materials is be- 
cause of their high density and high Z which is a de- 
sirable attribute for a scintillator material used in 7- 
ray detectors. [22 Ce can be readily incorporated as a 
dopant in many La and Y based scintillator materials, 
thus it is interesting to study the electronic structure 
of these materials doped with Ce. Another reason to 
study these compounds from a theoretical standpoint 
stems from some disagreement in the literature regarding 
La2Hf207:Ce. Borisevich et al. [23] reported Ce 3+ scin- 
tillation in La2Hf207:Ce powders. But, recent papers 
show no Ce 3+ luminescence for La2Hf2 07 when doped 
with Ce.[24][25] Ce 3+ luminescence in La 2 Hf 2 07:Ce is 
also counter- intuitive from simple empirical electroneg- 
ativity considerations [26 since electronegativities of La 
(1.1) and Ce (1.12) are smaller than Hf (1.3) on the Paul- 
ing scale. This would imply that Ce bd character states 
should lie above Hf 5d states. However, the electronega- 
tivity argument does not capture physical considerations 
such as local structure around the dopant, bonding, crys- 
tal field splittings, etc. In particular there are many cases 
of Ce doped La systems where the Ce 5d states lie above 
the La 5d states e.g. La203iCe and LaA103:Ce, even 
though the Ce has slightly higher electronegativity than 
La. [14] Ab initio calculations based on density functional 
theory have been successfully used by our group [27H29] 
as well as by others [30j [31] to study Ce 4/ and bd level lo- 
cations in known scintillator and non-scintillator materi- 
als. In this paper, we present results from first-principles 
bandstructure calculation of Ce doped RE2M2O7 class 
of compounds. The relative position of Ce 4/ and 5d 
with respect to the host valence and conduction band 
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edge, are computed to determine the possibility of Ce 3+ 
luminescence in these materials. 
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II. THEORY 

Scintillation in Ce doped materials corresponds to a ra- 
diative transition from the Ce excited state ([Xe]4/ 5c? 1 ), 
usually referred to as (Ce 3+ )*, to the ground state 
([Xe]4/ 1 5(i ). Our theoretical calculations for the pre- 
diction of candidate scintillator materials are based on 
the calculation of the Ce 5d and 4/ levels relative to the 
valence and conduction band of the host material. The 
Ce 4/ and bd levels must lie in the forbidden energy gap 
of the host material for Ce-activated luminescence to be 
observed in the material. 

The basic model for scintillation in a Ce doped mate- 
rial is that an incident gamma ray will produce a large 
number of electron-hole (e-h) pairs in the host material 
that then transfer to the Ce site. However, when the 
Ce 4/ level is in the valence band (Figure 1(a)) then a 
hole created by the 7-ray cannot be captured by the Ce 
and it may migrate to the top of valence band (VBM). 
This may result in fluorescence from the host but not 
Ce 3+ scintillation. On the other hand, if Ce 5d states 
are inside the conduction band then the electron created 
by the 7-ray will not diffuse to the Ce site since that is 
not energetically favorable (Figure 1(b)). Other possible 
scenarios for reduced or a lack of detectable Ce 3+ scintil- 
lation arise from loss of e-h pairs to competing processes 
in the host material [32, 33 (e.g., hole trapping, electron 
trapping, STEs) before transfer to the Ce site. 

Using first-principles calculations we can predict the 
Ce 4/- VBM energy gap in a number of host materials, 
and by studying the character of the lowest occupied con- 
duction band state we can qualitatively estimate if the 
Ce 5d states are below the conduction band. [27] We do 
not model other competing processes in the host which 
may reduce brightness. However, it should be noted that 
for Ce 3+ luminescence to be observed the Ce 4/ and bd 
states must "fit m"the bandgap of the host. This is a 
necessary (but not sufficient) condition for Ce 3+ lumi- 
nescence. 



III. COMPUTATIONAL PROCEDURE AND 
DETAILS OF CALCULATIONS 

Our calculation procedure follows the prescription de- 
scribed in Ref. [14]. We begin by taking atomic positions 
and crystal symmetry group information of the host com- 
pound from the Inorganic Crystal Structure Database 
(ICSD).[34j [35] The calculations are then performed in 
the following three stages. 

1. The host material is doped with a single Ce atom 
and the atomic positions are relaxed keeping the 
simulation cell volume and dimensions fixed. We 
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FIG. 1: Simple schematic model to illustr ate t wo 
possible scenarios for no Ce 3+ luminescence, (a) Ce 4/ 



below the top of host valence band |(b)| lowest Ce 5d 
inside the host conduction band. 



used a large supercell in order to minimize the in- 
teraction of the Ce 3+ dopant with its periodic im- 
ages. From computational considerations, a super- 
cell matrix containing 88 atoms was chosen for all 
calculations. By comparing with different cell sizes, 
this size of cell was found to give reasonably well 
converged results for the properties of interest. A 
single Ce 3+ ion replaces 1 in 16 ions of Y 3+ /La 3+ 
in the doped supercell. 

2. A ground state band structure calculation is per- 
formed to determine the position of the Ce 4/ level 
with respect to the VBM of the host material. The 
LDA+U approach [36] has been used to give a bet- 
ter description of the localized Ce 4/ state. 

3. A constrained LDA (excited state) calculation is 
then performed by setting the occupancy of Ce 
4/ states to zero and filling the next highest 
state. Subsequent analysis of the band decomposed 
charge density is used to qualitatively determine 
the position of Ce 5d relative to the bottom of the 
conduction band (CBM). 

All calculations were performed using the plane wave 
basis set code VASP.[37 39 Spin polarized calculations 
were carried out within the Perdew-Burke-Ernzerhof 
(PBE) parameterization to the exchange-correlation 
functional. [40] and the projected augmented wave (PAW) 
approach, [41 as implemented in VASP,[42] was used to 
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treat the valence electron-ion interaction. The plane 
wave basis set for the electronic wave functions was de- 
fined by a cutoff energy of 500 eV. Integration within the 
Brillouin zone was performed on a T centered 2x2x2 
grid of Appoints. The total energy convergence criterion 
was set to 10 -6 eV. Structural relaxations of atomic po- 
sitions were carried out until the residual forces on any 
atom were less than 0.01 eV/A. 

We have used the rotationally invariant method of Du- 
darev et al. [43] as implemented in the VASP code [44] for 
PBE+U calculations. We have used a value of U e ff = 
2.5 eV for the Ce 4/ states in these compounds where 
the Ce impurity is surrounded by nearest neighbor oxy- 
gen anions. This is based on our first-principles studies of 
Ce-doped compounds. [T4] The value of U e ff = 2.5 eV was 
found to be suitable for matching the Ce 4/-VBM gap to 
experiments for a variety of Ce doped compounds includ- 
ing oxides. DFT-PBE (or LDA) calculations incorrectly 
position the unoccupied La 4/ states at the bottom of 
conduction band as, for example, in LaBr3:Ce.[45 How- 
ever, La 4/ states lie higher in energy [46 so following 
the prescription of Ref. [47] we used a U e ff parameter for 
La 4/ states and move them higher in energy. An in- 
teresting aspect of RE 2 M 2 7 (RE = Y, La; M = Ti, Zr, 
Hf) compounds is that the transition metal ions have the 
same +4 valence, but with increasing principal quantum 
number the d-orbital becomes increasingly extended in 
character[48] from Ti (3d) to Zr (4d) to Hf (5d) com- 
pounds. Ti 3 c? states are most localized in nature and 
hence ill described by LDA (or GGA) because of the as- 
sociated self-interaction error. [49 Our PBE+U calcula- 
tions include a +U correction for Ti 3d states with U e ff 
= 5.5 eV.[50] 

A constrained LDA calculation is performed at the T 
point by manually setting the occupation of the Ce 4/ 
states to zero and filling the first excited d character 
state. From the analysis of the band decomposed charge 
density we determine the localization of the excited state 
on the Ce site by measuring the percentage of this state 
surrounding the Ce atom. This is done by constructing a 
Voronoi type polyhedron centered on each ion and inte- 
grating the portion of charge density within this volume. 



IV. RESULTS AND DISCUSSION 

A. Ground State 

The La compounds studied in this work mostly crystal- 
lize in the cubic pyrochlore structure (space group Fd3m) 
except La2Ti207, which exists in a monoclinic layered 
perovskite structure (space group PI 121) at room tem- 
perature. The pyrochlore structure[51 (A2B2O7) is a 
super structure derivative of the simple fluorite structure 
(AO2 = A4O8), where the A and B cations are ordered 
along the <110> direction and one-eigths of the anions 
are absent. The anion vacancy resides in the tetrahedral 
intersticial site between adjacent B-site cations which re- 



duces the coordination of cation B from 8 to 6. The 
pyrochlores are generally formed by combining a triva- 
lent A cation and a tetravalent B cation. The monoclinic 
structure is a variant with same stoichiometry as the py- 
rochlore structure but it is preferred when cation A is 
larger than cation B. 
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FIG. 2: Total and atom projected partial density of 
states plot for La2Ti207:Ce from PBE+U spin polarized 
calculations (U e ff = 2.5eV for Ce 4/ states, U e ff = 
10.36eV for La 4/ states and U eff = 5.5eV for Ti 3d 
states). Spin up states are shown in solid lines and spin 
down with dashed lines. Ti 3d states are at the bottom 
of conduction band while Ce 5d states are higher in 
energy. Oxygen 2p forms the top of valence band. 

Figure [2] shows the total and atom projected partial 
density of states (DOS) for monoclinic La2Ti 2 07:Ce. 
The Fermi level is set at eV which is also the po- 
sition of the occupied Ce 4/ state. Oxygen 2p states 
mainly form the top of valence band and the host VBM 
is separated from the occupied Ce 4/ (i.e., Ce 4/- VBM 
gap) by approximately 2.5 eV. We can see clearly from 
the ground state DOS plots that the conduction band 
minimum (CBM) has a predominantly 3d character at- 
tributed to Ti. Ce (and La) 5d states are about 2-3 eV 
higher in energy. Therefore, it is highly unlikely that an 
ionized electron will localize on the Ce site. An impor- 
tant point to note here is that the nature of the CBM is 
markedly different from known oxide Ce scintillators such 
as LaB30e:Ce and YA10s:Ce, where La or Y d charac- 
ter states are generally the lowest lying host conduction 
band states. 

Amongst all the compounds studied in this work, only 
La2Ti207 has more than one inequivalent trivalent sites 
(La 3+ ) for Ce 3+ substitution. These La sites differ in 
their oxygen coordination and bond lengths. We found 
that the Ce 4/- VBM gap varies between 2 eV to 2.9 eV 
for the different Ce 3+ substitution sites (Figure |3|. How- 
ever, since the Ce 5d is situated well above CBM, even 



though the crystal field effects change due to change in 
local environment, Ti 3d states are still at the bottom of 
conduction band. Figure [2] shows partial atom projected 
DOS plots when Ce is substituted on La site III. 
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FIG. 3: Total density of states plot for La2Ti207:Ce 
for different Ce substitution sites from PBE+U spin 
polarized calculations. Spin up states are shown in solid 
lines (brown) and spin down with dashed lines 
(magenta). 

Figure [4] shows the partial DOS plots for I^Z^OyiCe 
and La2Hf2 07:Ce. Oxygen 2p states predominantly form 
the top of valence band. Because the host compounds 
have same cubic crystal structure, the Ce dopant envi- 
ronment is similar (surrounded by same number of near- 
est neighbor oxygen atoms) in the doped systems. The 
4/-VBM energy gap are therefore very similar for these 
systems. As we move from Zr to Hf compound, we see 
that the host conduction band states are progressively 
an admixture of La and transition metal d states. This 
systematic trend of overlapping of transition metal d 
states with La 3+ d states has been analogously observed 
for SrMOs (M = Ti, Zr, Hf) oxides from spectroscopic 
measurements. [52] 

A common feature of these ground state plots is the 
positioning of unoccupied Ce 4/ states at the bottom of 
the conduction band. However, in a situation of hole cap- 
ture by the Ce ion, Ce 4/ states will shift and move lower 
in energy. We have seen this movement of unoccupied Ce 
4/ states (sometimes by about 2-3eV) from excited state 
calculations. [14] This is due to the localized nature of the 
Ce 4/ states in these compounds. Considering the states 
above the Fermi level, for Ti and Zr compounds, Ce 5d 
character states are located clearly above the bottom of 
conduction band. However, Ce 5d states show mixing 
with the host d character bands at the bottom of con- 
duction band for La2Hf2 7 :Ce. This is noteworthy since 
the nature of conduction band states in this compound 
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FIG. 4: Atom proje cted partial densit y of s tates (DOS) 

plots for Ce doped (a) La2Zr207 and (b) La2Hf207. 
Fermi level is set at OeV. Spin up states are shown with 
solid lines while spin down with dashed lines. Oxygen 
2p states form the top of valence band. 



are very similar (5d states of host ions La and Hf and Ce 
dopant) and we expect the DFT-PBE error in relative 
positioning of these states to cancel to a large extent. It 
is, therefore, important to check the localization of the 
excited state on the Ce site for La2Hf20 7 :Ce to estimate 
the possibility of a Ce 5d-4/ transition. 

We also studied Ce doped compounds of Y i.e., 
Y 2 M 2 7 :Ce (M = Ti, Zr, Hf). Y 2 Ti 2 7 crystallizes 
in cubic pyrochlore structure however, Y2Zr2 7 and 
Y2Hf2 7 prefer the closely related disordered fluorite 
structure. The disordered fluorite structure exhibits the 
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Fm3m space group and has 50% A and B cation occu- 
pancy on the cation sublattice and a 7/8 occupancy on 
the 8c oxygen sublattice. This structure has the same sto- 
ichiometry as the pyrochlore structure, but is preferred 
when cation B is larger than cation A. Disordered fluorite 
structures of Y2Zr207 and Y2Hf2 07 are not listed in the 
ICSD so our calculations were limited to the related, but 
unstable, pyrochlore phase for these compounds. Fig- 
ure [5] shows the total and atom projected partial density 
of states plot for Y2Ti207:Ce. Ground state DOS plots 
for this material show features similar to I^T^C^Ce 
with Ti 3<i states forming the bottom of the conduc- 
tion band and O 2p predominately at the top of valence 
band. Ce 5d character states are situated much higher 
in the conduction band, which indicates that this ma- 
terial will not show Ce 3+ luminescence. Similar to La 
compounds, systematic shift in the position of transition 
metal d states was observed in Y2M20r:Ce with Zr Ad 
and Hf hd showing progressive overlapping with Y Ad 
states. 
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FIG. 5: Total and atom projected partial density of 
states plot for Y 2 Ti 2 07:Ce from PBE+U spin polarized 
calculations (U e ff = 2.5eV for Ce Af states and U e fr = 
5.5eV for Ti 3d states). Fermi level is set at 0. Ti 3d 
states form the bottom of conduction band while 
Oxygen 2p form the top of valence band. 



B. Excited State 

For most Ce-doped scintillator candidate materials, Ce 
5d states hybridize with the host conduction band states 
like La 5d or Y Ad states. A constrained LDA calcu- 
lation gives us a qualitative measure of the position of 
the Ce 5d state relative to the CBM and from our ex- 
perience these calculations give consistent results for Ce 
doped scintillators and non-scintillators.[l4] However, for 



the class of the systems studied in this work, because the 
bottom of conduction band is mostly composed of the 
transition metal d character states, it is expected that 
the localization would be mostly around the transition 
metal atoms. 

Figure [6] shows the charge density isosurface plots of 
the first d character conduction band state for the com- 
pounds studied in this work. Localization on the Ce site 
would indicate possibility of Ce 3+ luminescence whereas 
a delocalized state indicates no possibility of Ce 3+ lumi- 
nescence. We can clearly see that in all compounds stud- 
ied there is no localization of the excited electron state on 
the Ce site. Therefore, our calculations indicate that the 
Ce 5d states are situated much higher in energy in the 
conduction band and, as a consequence, ionized electron 
will not localize on the Ce site. Thus, Ce luminescence 
in these host materials is unlikely. 

It is important to note that Ce 3+ activated AHfOs 
(A = Ba, Sr, or Ca) ceramic scintillators [53| [54] have 
been reported in literature. These systems are however, 
quite different from the materials studied in this work. 
Foremost, there are no trivalent sites available for natural 
substitution of Ce 3+ (Ba, Sr and Ca are divalent while 
Hf, Zr and Ti are tetravalent) unlike La or Y compounds. 
Also, understanding the nature of charge compensation 
if, for example, Ce 3+ substitutes a Ba 2+ , requires more 
work such as shown in reference. [55] Consequently, these 
hafnates are not subject of our present study. 

During the course of our investigation we found 
more examples of La and Y compounds with Ti, Zr 
or Hf as constituents which can be easily identified 
as unsuitable candidates for Ce 3+ luminescence from 
first-principles calculations. Ba6Y 2 Ti 4 0i7, LaLiTi04, 
LaNaTi0 4 , YNaTi0 4 , LaZrF 7 , Y 2 HfS 5 are examples of 
predicted non-scintillating Ce-doped host materials. To 
our knowledge there are no Ce 3+ luminescent La or Y 
host materials (phosphors or scintillators) with Ti 4+ and 
Zr 4+ as constituent ions, which agrees with our calcu- 
lations. As far as La2Hf207:Ce is concerned, our cal- 
culations agree with the observation of Ji et al [24] and 
Cherepy et al. [25] in that we believe Ce luminescence in 
this material is unlikely; although further experimental 
data on this material will be helpful. This material is the 
subject of experimental investigation by our group and 
detailed results will be published in the future. We would 
like to also mention that a recent experimental paper on 
Lu 4 Hf30i2 powders doped with different rare-earth ions 
did not find any evidence of Ce activation in this host. [56] 
Havlak et al.[56 think that this is probably due to the 
positioning of lowest Ce 5d inside the host conduction 
band. This result is significant since Lu 4 HfsOi2 crystal- 
lizes in the (5-phase (space group R-3m) which is closely 
related to the fluorite structure. The pyrochlore and 5- 
phase are oxygen deficient fluorite structure derivatives. 
Further information on these structures can be found, for 
example, in Ref. [57] . 

Luminescence of Ti 4+ in La2Hf207 host has been stud- 
ied for application in X-ray imaging. [20] Our calculations 



6 




(d) (e) (f) 

0Ce 9 La $Ti Zr Hf $ C 



FIG. 6: Excited state plots for RE2M207iCe compounds (RE = Y, L a; M = Ti, Zr , Hf ). Plots s how charge de nsity 
isosurface s of the lowest exc ited d state at 50% iso-surface threshold. |(a)| Y2Ti2Q7; [(b)] Y2Zr2p7; |(c)| Y2Hf2p7; |(d)| 
La2Ti2Q7; |(e)| La2Zr2Q7; |(f)| La2Hf2Q7. No localization on the Ce site indicates no Ce scintillation. Coloring scheme 

is indicated in the figure. 



for La2Hf207iTi clearly show the presence of Ti 3d states 
well below the bottom of conduction band composed of 
an admixture of La and Hf 5d states. This indicates that 
a charge transfer transition between the O 2p states and 
Ti 3d states is possible in agreement with measurements. 

Figure [7] further examines the three lowest conduction 
band states of La2Hf 2 07iCe from excited state calcula- 
tions. From our calculations we scan up the higher ex- 
cited states and search for states which may have pre- 
dominately Ce 5d character. We find that, although the 
lowest excited state has a mainly Hf bd character, the 
next highest states ((z+1) and (i+2)) are mainly local- 
ized around the Ce site. This situation is analogous to 
the results presented by van der Kolk et a/. [58 from ex- 
perimental studies of LaAlOaiCe where they place low- 
est Ce 5d states few tenths of an eV above the bottom 
of conduction band. Our calculations ignore the effects 
of spin-orbit coupling on the splitting of 5d character 
states of the conduction band. A recent paper by Gra- 
cia et al.[59 presented results from cluster calculations 
of Ce-doped Y 3 Al 5 0i 2 (YAG). Ce 3+ substitutes for Y 3 + 
ion in an eightfold coordination of oxygens in YAG. Their 



results show that, upon including the effect of spin-orbit 
coupling the Ce 5d states moved by about 0.1 eV higher 
in energy. Meanwhile, Garcia et al.[60 find that with the 
inclusion of spin-orbit interactions, Hf 5d states shift by 
about 0.05 eV in Hf0 2 using FLAPW calculations. This 
suggests that including the effects of spin-orbit coupling 
will not change the relative position of the lowest Ce 5d 
state with respect to the bottom of the conduction band 
(Hf 5d character). 



Results presented so far have demonstrated one possi- 
ble situation for no Ce luminescence wherein the lowest 
Ce bd is above the CBM (Figure 1(b)). This is also the 
reason why compounds such as LaAlOaiCe, La203iCe, 
etc, do not show Ce 3+ luminescence, and our calculations 
show that in agreement with experimental results. [14] In 
our studies of more than hundred Ce doped compounds 
so far we have not seen evidence of Ce 4/ situated in 
the host valence band of candidate scintillator materials. 
This is consistent with available experimental data on Ce 
doped compounds. 
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FIG. 7: Charge density iso-surface plots of the three lowe st ex cited states of La2Hf20r:Ce at the gamma point, (a) 
lowest excited s tate (state i) with mainly Hf 5d character; |(b)| next highest excited state (state z+1) with mostly Ce 
5d character; |(c)| next higher excited state (state z+2) with predominately Ce 5d character. State i is separated 
from state (z+1) by 40 meV while states (z+1) and (z+2) are separated by 6 meV. 



V. CONCLUSION 

We have presented results of first-principles density 
functional theory calculations for Ce doped RE2M2O7 
compounds (RE = Y, La; M = Ti, Zr, Hf). Our cal- 
culations indicate that the lowest Ce 5d state lies above 
the bottom of conduction band in these materials. This 
would prevent Ce luminescence in these materials. Our 
results agree with available experimental data. 

The predictive power of any pre-screen selection cri- 
teria is critically dependent on the successful prediction 
of negative candidates. This paper presents results of a 
class of materials which are found to be non-scintillators. 
These results indicate that our systematic calculation 



procedure can successfully screen non-scintillators and 
therefore, be applied to effectively identify candidate ma- 
terials for Ce 3+ luminescence in inorganic compounds. 
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